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The motion of a charged particle in a 2-dimensional neutral 
sheet with linear magnetic field variation is analyzed by Hamil- 
tonian methods. From the basic Hamiltonian formulation, results 
of Sonnerup (J. Geophys. Res. 76 , 8211, '7l) are recovered and 
given intuitive interpretation. The transformed Hamiltonian is 
used to derive the correct frequency of oscillation and serves 
as the basis of analytical treatment of perturbed versions of the 
motions e.g. with electric field added. For 2-dimensional fields 
with slightly differing configurations - in particulars with added 
small field component orthogonal to the sheet, with an X-type null 
or with a chain of null points -two alternative methods are deve- 
loped, reducing the problem either to motion in a 2-dimensional 
potential or to a pair of coupled oscillators. 



Charged Particle Motion In the Vicinity of a Neutral Plane 


D. Stem 1.19.72 


The Hamiltonian 

Let a magnetic field B be given, directed along the cartesian y-axis 
with the field intensity proportional to x and vanishing on the plane x - 0 

B = (Bq/l) x £ (1) 

This somewhat resembles the field in the tail of the magnetosphere: the plane 
x ss 0 corresponds to the neutral plane and L is a scale length, giving the 
distance from the neutral plane at which the field intensity reaches some 
standard value B Q (The notation so far corresponds to that of E.N. Parker, 

Phys. Rev. 107. 924 (1957), who included motion in such a field as an example 
of non-adiabatic motion at the end of his article on the guiding-center plasma; 
the cartesian axes, however, correspond to the choice made in "Adiabatic 
Particle Orbits in a Magnetic Null Sheet" by B.U.fl. Sonne rup, JGR 8211, '7l) 

By introducing Euler potentials ( , f ) one easily derives a vector 
potential A of the form : 

B - (- x B q /L ) ( 7 x X V z ) = - (Bq/ 2 L) ( Vx 2 X Vz ) 

= V * (-x 2 B 0 /2 L)i (2) 

The Hamiltonian is thus 

H “ k{ P x 2 + p y 2 + ( P z + («B ( /2 Lc) x 2 ) 2 j- (3) 

Two things at once stand out. First of all, the motion separates into an 
unconstrained motion along the field's direction (y-axis) and a motion 
orthogonal to it. If we abbreviate 
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eB ( y / 2cL =0" ( 4 ) 

we get for the latter motion 

Hl = (l/2m) | p x 2 + (p z + Q*x 2 ) J » = const. ^ 

Secondly, since z is absent, p z is a constant of the motion which we 
will denote by P g . In terms of the particle's velocity 

. 2 < \ 

P z = mz + (j"x = const. (6) 

It is also convenient to define a constant Pj_ equal to the momentum component 
orthogonal to B 

p i = 2 m B x (7) 


Qualitative Aspects 


Because of th© constancy of P g , the motion produced by the Hamiltonian ( 5 ) 
may be viewed as a one-dimensional motion in a potential 

V - (l/2m) (P a + <Tx 2 ) 2 (8) 

Two possibilities exist, depending on whether P z and CT have identical or 
opposite signs. In the first case we may define, a suitable constant Xq and 
write (note that V is in any case positive-definite) 


V = ( <I 2 /2m) ( r 2 «■ x 2 ) 

0 (9) 

In that case it is found that V has a single 
minimum at x = 0 , and its value there is 

Xq 2 /2e . The potential well is then parabola- 
shaped (but given by a 4th degree curve) as 
shown in the sketch: 




In this case the particle oscillates 
from one side of the well to the other. 
With each oscillation, it crosses the 
neutral plane x = 0 : the motion is 

therefore as drawn in the sketch on the 
right. Following Sonnerup we will call 
this a meandering mode of motion — 
specifically, meandering without looping; 
itself at any point. 

If P z and <T have opposite signs, i 
so that the potential assumes the form 



suitable constant Xq may be defined 


V = ( <T 2 / 2 m) ( x 2 - Xq 2 ) (10) 

7 then assumes its smallest possible value, namely zero, at x = * x^ . 

If one examines V/^x one confirms that these are true extrema, and 
that there also exists an extremum 
at x = 0 9 which clearly must be 
a relative maximum. The graph of V 
therefore has the general appearance of the 
curve drawn on the right. 

Two possibilities exist: if the 
particle is trapped in one of the 
"side pockets", it will oscillate 
back and forth without crossing the 
neutral sheet x = 0 . The motion then has the general nature of guiding-center 
motion, especially for orbits far away from the neutral sheet. For orbits that 
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come close to the neutral plane 
(i.e* those oscillating rather 
high up in the "pocket") there 
will exist a marked asymmetry: both 
cases are shown in the drawing on 
the rights ^ 



If the particle clears the central peak, on the other hand, it will again 
oscillate across the plane x = 0 and again the result will be a meandering 
mode of motion. This mode, however, 
will have orbits crossing themselves 
("meandering with looping"). To show 
that this indeed occurs, we us© the 
Hamiltonian ( 5 ) to derive z . We 


a =-( <r / m) ( x t * 0 ) 



inner 


the’ sign being positive if and 

are of the same sign and negative if they are not. In the first case, z is 
always of the sane sign, i.e. the particle steadily advances in one direction 
on the E-axis, as in meandering without looping. In the second case, z can 
assume either sign — it changes sign as the particle passes x = ± Xq — and 
this implies looping. 


As Sonnerup has shown, there exists a continuum of solutions ranging from 
regular guiding~center modes to meandering modes, in the field given here. 

One may look at the situation from the following angle. Far from the neutral 
sheet the motion is of the guiding-center type and the particle's drift velocity 
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has the order of magnitude 

v drift constant • ( pf |vb|/c*>B ) (12) 

where 0) is the angular velocity of gyration and all quantities are 
evaluated at the guiding center. In our particular example, the enumerator 
is constant ( V B has only one component, equal to Bq/L ) while each 
term in the denominator is proportional to B . Thus as the guiding center 
approaches the neutral plane, the above expression diverges as B . 

Physically, the formula for the drift velocity is expected to break 
down somewhat before this situation is reached, because a particle cannot 
drift at a velocity larger than some fraction of P^/m . This indeed 
happens when the guiding-center mode of motion evolves into the meandering 
mode, in which the z-component of the velocity indeed is some appreciable 
fraction of P|./m . 

Analytical Treatment 

Let us follow the Hamilton-J acobi method for generating a canonical 
transformation 

( z, p 2 , x, ) -*• ( Pj., Qg, P 2 ) (13) 

such that the new momenta are constants and the new coordinates are linear 
in tim©o Since p z already is a constant we may write the generating 
function, tentatively, as 

W - z P x + V^, P 2 , x) (14) 

The Hamil ton-J acobi equation then gives (P^ = P z , obviously) 

(^/Dx) 2 + (P 1+ C rx 2 ) 2 * Pi ( 15 ) 
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This contains Pj_ wh ich depends directly on the constant of energy E x . 
In general, P x will also be some function of (P^, P 2 ) ; this functional 
relationship may be specified more or less arbitrarily and every choice 
leads to a different definition of P 2 • Any of these definitions leads to 
a solution of the Hamilton-J acobi equation; however, if we plan to deal with 
problems differing from the given one by some small perturbation, it is 
advantageous to make a special choice, namely, to define P 2 as the action 
variable J , evaluated as an integral over one period of the system: 

J - K 4 * - flf 4 * 


1 


P? 


(Pi 


+ <T x 2 ) 2 


dx 


(17) 


This choice arises naturally if on© demands that its canonical conjugate 
■ — denoted in this case by XI and called the angle variable — increases 
by unity each period: 

1 ■ ,f d!1 “It ! 41 ■ f ih d * 

- f & t: 41 * (18) 

In cur case this leads to a general type of solution 


<f P s to « J + P( P 1# <T ) (19) 

but it will b© presently seen that it is advantageous to set the added 
arbitrary function equal to zero. Intuitively, this choice of means 
that when the Hamiltonian is expressed in terms of the new variables, XI 
will always appear in functions that are periodic in XI with period unity, 
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that is, sin(27T£l.), cos(271£i) and their higher harmonics, and functions 
which can (in principle) be expanded in fourier series in such functions. 

Mathematically, the result may be viewed as follows. When all dependence 

on is of the type described above, then any function U of the variables 

can be separated into the sum of two terms: the "average" of U (over O. ) 

1 

<”>= 5 U da s J U dfl (20) 

0 

which does not depend on fY at all, and the "purely periodic" part of U 

U per = U " Ul) 

which has zero average and therefore represents a purely oscillating 
quantity. 

The significance of the choice of J as in (17), with the additional 
function allowed by (19) set equal to zero, is that when it is used, all 
partial derivatives of are purely periodic, except for ^W^/t> J 
which increases by unity each period. As an example, consider the derivative 
by P^ % for this to be purely periodic, one requires 

" 0 = - f ’x 1 * 

( 22 ) 

" ^ ■ 0 

The last equality follows from J and being independent (new) variables. 
The fact that 'D/'DP^ may be taken out before the integral follows from 
the fact that the integration is over a closed loop: it could be broken up 
into two or more sections, and the limits of integration in each section 
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would, then be spelled out. These limits may well depend on P^ , but it is easy 
to see that the contribution of each limit appears then in two terms that 
cancel each other, an upper limit term and a lower limit term* 


Some consequences of this may be pointed out. First of all, will 
differ from z only by purely periodic terms, for 




w _ j) w i 

^ p i = 2 + Vi 


(23) 


and the second term is purely periodic. 


Next, suppose 0~ is a slow function of the time t . Then a term 
is added to the Hamiltonian, in the new variables, equal to 


2v dcr ^> w i 

^>t d t 


(24) 


If we are to handle the motion by perturbation methods', the addition to the 

Hamiltonian must be small, with smallness defined by some small parameter £ . 

time , the l ight side of j 

If the variation of 0” (t) is "slow", then the^derivative on r (24) is indeed 
of order £ , but it is equally important to note that the factor that multiplies 
it is "purely periodic", for if it were otherwise, the term would gradually 
grow until it is no longer "small". 

A similar situation occurs when O’ is a slow function of z . Now p 

z 

is no longer a constant of the motion, but we can apply the transformation 
(16) anyway and get 

p z = 'Dw/'frz m p x + (dcr/dz) ( "1 Vj/lXT) 

(25) 

=» P^ + (purely periodic function of order £. ) 


This turns out to be necessary for a successful perturbation solution here 
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Calculation f or t h e Zero-order Case 


In the present case the integral (17) — which essentially gives the 

complete 

relationship between P, and J — may be expressed in terms of. elliptic 

**• A 

integrals of the first and second kind 

3T/2 


K(k) 


du 


i VT7 


2 2 
k sin u 


\ ( 1 - J * 2 * 64>= 4 + (*> 


E(k) 


*/2 

J if 


2 2 

1 - k sin u du a 


I - 4 k2 - 64 k4 — > (J7) 


In order to reduce the integral (17), the square root must be first transferred 
to the denominator. Two ways of achieving this are possible • We may integrate 
by parts (note that there is no contribution from the limits since integration 
is around a closed path) 

dv a - (j)" 

Choosing 


v du 


21 * 


v 3 r 


U 3 j jJ - (P x + X 2 ) 2 } 


this gives 


^ 2 P 1 <Tx 2 + 2 O' 2 x 4 
u 


dx 


Alternatively, we may write 


J 3 I, 


£ ax = 4> 


r 2 2 2 2 4 

pf - P, - 2 P, (Tx - <r x 4 


u 


J 


( 28 ) 


(29) 


u 


dx (29) 
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4 _ 

To eliminate x from the enumerator, we express J as 



where we define 

yW 2 = ( Px - P x )/c r (31) 

This will always be positive, by the following arguments. From the form of the 

Hamiltonian, nothing is changed if we simultaneously reverse the signs of both 

CT and ; we may therefore assume that CT is positive* The form of the 

2 2 

Hamiltonian also shows that ¥ x exceeds P^ , and since Pj_ is defined as 
a positive quantity, it follows that the expression on the right side of (3l) 
is also positive. 

One now substitutes 



* Note: z used in (32) and in the equations that follow is an auxiliary variable 
not related to the z-coordinate. 
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Let be eliminated and replaced by k , defined by 

k 2 = i {l - (Pj/Pj} (34) 

ioOa 

? x = P x ( 1 - 2 k 2 ) (35) 

The denominator in (33) then becomes 

+ Pj^ + (^ - P 1 ) cos 2 z = (2P X - (P A - P x ) sin 2 z 

» ^2 ~P^ |l - k 2 sin 2 a]'^ (36) 

The mmierator is 



Substituting everything and noting that the range of integration of y is 
4 times that of K(k) or E(k) ( 0 to 2JT instead of 0 to lt/2 ), we get 
(Sonnerup 9 s notation): 



Prom general perturbation theory it follows that J is conserved adiabatically 
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if O' varies slowly with z and/or time. As an example demonstrating the 
intuitive significance of such conservation, consider a particle in the 
meandering mode with small P^ • Then 

k 2 % £ f x (k) = JT/4 = constant (39) 

J 2 = Px/CT (40) 

Now by eq. (4) O' is proportional to Bq/l, which is the only non-zero 
component of VB . If (T grows slowly in time, Px also grows slowly, at 
a rate that leaves J approximately constant. 

If the meandering motion carries the particle to regions of stronger 
gradient (i.e. with 0*(z) ), J is ©till an adiabatic invariant® 

However, since energy is conserved and the Hamiltonian (it can be shown) 
still separates as before, Pj. cannot vary and the conservation of J 
is maintained by the variation of k , through changes in P^ . Equation 

(40) may not be used in this case. 

The New Hamiltonian 

Since k is a function of and P^ , equation (38) defines an implicit 
relationship 

p j. - ? X (J. (41) 

When this is substituted in (7), it gives the new Hamiltonian (which without time 
dependence equals the old one). The trouble is that inverting (38) to the form 

(41) is not easys for instance, we could use the series expansions (26) and (27), 


but this tends to be cumbersome • 
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There exists a nay around this difficulty, namely by redefining P^, the 
new variable introduced by the canonical transformation generated by W . 
Instead of defining it (as before) equal to p z , we define 

P x = k » iJ i [l - (P/Pjl (42) 


Equation (16) is now replaced by 



and (42) follows immediately from 

p » OW/njz (44) 

z 

Also, since 

P x = '~bV 1 /'c) x (45) 

the derivation of J is unchanged, giving instead of (38) 

J » (16/3) P 3 / 2 (2 <T)“ * f^) (46) 

from which 

H - (l/2m) (3/16) 4 / 3 (2(T) 2 / 3 J 4//? {^(P].)}^ 3 (47) 

Since jQ. grows by unity each period, the frequency V of the oscillation 
is at once obtained from 


y - 


(2/3m) Pf/J 


(48) 
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This choice of will be termed the second choice . It is useful for 
perturbed motions where the perturbation consists of some additional terms 
in the Hamiltonian (examples of this will be described later) but not where 
<J" is a slow function of position or time. Suppose that (T » CT(t) ; then 
the new Hamiltonian acquires two added terms 


21 


z (l - 2 ? x 2 ) 


^Pj . 

r b(T 


d<T 

dt 


+ 


V<T 


d<r 

dt 


(49) 


Both terms contain a "small" factor dcr/dt . In the second term, the 
expression multiplying this factor is purely periodic, by arguments cited 
earlier, but in the first term we get an expression proportional to z , 
and this is liable to grow as the motion progresses. In such a case, there- 
fore, we must use the first choice of P., equal to p and invert (38) 
to (4l) by other means. 


It is also possible for to depend on x , which results in a different 
shape of the oscillation. The problem then is rather similar, but the 
integrals one obtains cannot be reduced to elliptic form and must be handled 
by expansion, numerically or by other means. 


Slow Dependence of CT on z 

As an example, assume that <T = d(a) , with dcr/dz small compared, say , 
to the ratio between CT and the distance between two consecutive crossings 
of the neutral plane (in the meandering mode). The Hamiltonian again separates 
into two parts, and the part describing the motion orthogonal to the field 
(we will drop the subscript here) is 

H = (l/2m) £ p 2 + ( P 8 + <TU) x 2 ) 2 } (50) 
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Obviously p z is no longer conserved, but we apply anyway the transformation 
generated by (16). By (44) and (45) » with the "first choice" of P^ and with P x 
defined as in (58) 

P x = ^W/'Dx - P 2 - (P x + <Tx 2 ) 2 (51) 

p * 'Dw/'&z = P, + ( / T>w/'d(T) d<T/dz (52) 

Z X 

Substituting this in (50) gives as the new Hamiltonian 

H (l/2m)[pf - (P 1+ (Tx 2 ) 2 + (P x 4 + <rx 2 ) 2 1 

= (l/2m)|pf + 2 ST (Pi+cTx 2 ) + 0(£ 2 )]' (53) 

The zero-order part of this is a function implici tely defined 

through (38) 

P l /2m = H (0) (J, P 1 , <T( z)) (54) 

It is not permissible, however, to use z in the new Hamiltonian, and we 
therefor® transform and expand, using (23) 

o w 

cr(x) = cr^) - + o(e 2 ) (55) 


Because (f(z) differs from <r(Q 1 ) only by small terms we may replace 

in the first-order part of (53) d<r/dz by d <T /dQ^ . In the same 

2 

expression we also must replace x by a function of the new variables. 



Now the angle variable ^ s > with the same definition of z 
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Hence (e.g. eqs„ (133), ( 135) in section 112 of Franklin, Methods of Advanced 



The new Hamiltonian thus is 
H = H ( °\j, P r ct(q 1 )) 


^H^ dcr 2S(i 
2)(T d^ + 


1 2LS Ml 

m 'D CF 



f XlVscrp^ l \ 

V* pf/3 j J ) 


+ o(e 2 ) 


(58) 


The 0(6-) tern is purely periodic with zero average, since periodic factors 
only appear in it in odd powers, so J is indeed adiabatically conserved 

J a constant + 0( £ 2 ) (59) 

One could go on here and derive higher-order corrections, but there seems to 
exist little practical need for those. 
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Adding a Constant Electric Field Along z 


Suppose that we add a constant electric field in the z-direction 

B = - (60) 

4> ** - Eq z (6l) 

Then 

H = (l/2m) £ P x 2 + (p z + <T x 2 ) 2 } - 8 E 0 2 (62) 


If the electric field is relatively small, we can regard it as a perturbation 
superimposed on the type of motion previously analyzed. In beginning the 
calculation we disregard the perturbing term and apply the canonical transfor- 
mation (43), where (P^, J) are the new canonical momenta, P x is defined by 
(38) and the "second choice" (42) is used in defining P^ . The main part of 
the Hamiltonian transforms as before and we get 

H * pj /2m - e B 0 z (63) 

To express z one must use (23) 


Ql = 



22l 

D ?1 + ^ 


(64) 


where we have introduced an auxiliary function equal to 


P Z (J, P r cr ) » P x (l - 2P X 2 ) (65) 

Thus 


H * ( 3 / 16) 4 / 3 (2 O ') 2 / 3 J 4 / 3 [ f ^ Pj )] "^ 3 


e E, 


^7 


'DP. 



^ J 


( 66 ) 
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The Hamiltonian (66) contains 3 terms, only the last of which (the one 
proportional to ''DW^/'DP^ ) involves XL . If only that term is considered 
as a perturbation term, canonical transformations may be devised that push the 
dependence on XL to an arbitrarily high order* There remains a one-dimensional 
system in (P^, Q^) which is readily solved, at least in principle. 

Small Added Magnetic Field in x-direction 

Suppose that a small magnetic component 

£ = b 9 (67) 

is added, so that instead of a neutral sheet we have a sheet of very low field 

intensity and field lines become parabola-shaped 



Several approaches to this problem exist and will be described in what follows, 
(l) Intuitive Treatment 

It is easy to make a crude guess about what might happen here. No additional 
force exists in the x-direction, so the oscillatory motion orthogonal to the 
plane x =» 0 should not be affected. If this motion is averaged out by vise of 
action-angle variables (J, fl) one is left with an average motion in that plane, 
presumably one affected only by b . We could then get a circular motion, 
resembling the guiding-center motion due to b, alone : 
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Sonnerup , who treated this case numerically, indeed obtained certain 
solutions which resemble this form, though for other solutions a marked 
y - z asymmetry is evident (see Figure 4 in his paper). We may view 
the situation given in the sketch above as a motion with the visual two 
lowest invariants, but with the ordering of their characteristic time scales 
reversed: the oscillation about the orbit's point of minimum field inten- 
sity is rapid, while the guiding center gyration is slow. 

(2) Treatment as Coupled Oscillators 

Analytically, we have 

b = b: = b(VyxVz) 

= V x ( by z ) (68) 

The y-component of the motion can no longer be separated and the Hamiltonian 
thus is 

H = (l/2m) |p x 2 ♦ p y 2 + (p a + CToc 2 - (eb/c) y ) 2 | (69) 

On the other hand, p z is still a constant of the motion 

p = P 

*z z. 


constant ( 70 ) 
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Introducing new notation 

<Ti = O' <r 2 = eb/c 

we can rearrange H to the following form 


H 



(P. 


2x2 

x ) 


} * tf 


(P„ 



- P g 2 /2m - ( <T X <T^ /m) x 2 y (7l) 


This can be viewed as consisting of four parts. The first part is an 
enharmonic oscillator in the x-direction, of the type encountered in the 
treatment of meandering motion; the second term represents a harmonic oscillator 
in the y-direction j the third term is constant and does not affect the 
motion* while the last term represents an interaction between the two 
oscillators. It is proportional to b and therefore relatively small. 

To treat the motion we transform the Hamiltonian to the action-angle 
variables of the two oscillators and then treat the interaction term by 
classical perturbation theory. The problem resembles & wide class of problems 
in celestial mechanics and may involve resonance problems, since it involves 
the interaction of two oscillating systems with different frequencies. In 
particular, the problem of two harmonic oscillators coupled by a nonlinear 
term similar to the on© in (7l) has been extensively studied by Contopoulos 
(Astrophys. J. 153. 83, *68 ; Astronomical J. 73. 86, '68 ; Astronomical J. 

25, 96, 108, '70 ; Astronomical J. 21, 147, '7l). 
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( 3 ) Treatment as Motion in a Potential 

Because P z is constant, the Hamiltonian 

H » (1/2*) { p x 2 + P y 2 + (P z + <r x x 2 - <r 2 y) 2 } (72) 

may be viewed as describing motion in a potential 

V = (l/2m) (P z + (T^ 2 - (T 2 y ) 2 

y* 2 

= (cr 2 2 /2m) {(^/ay* 2 - (y-y Q )} (73) 

The potential V is never negative and from its second form in (73) it 
follows that for every value of V there exist two parabolic equipo tentials 
of the form 


y = y Q + ( ( J' 1 /cr 2 ) X 2 t m v/<r 2 (74) 

These are equally spaced in relation to the equi potential 

y = y 0 + ( cry cry x 2 (75) 

which corresponds to the lowest level V = 0 „ Some typical contours are 
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The topography of V is that of a parabolic channel which (when viewed in 

normal cross-section) narrows down with increasing y . For particles 

which move down the channel for* considerable distance, the usual two 

particles^ 

adiabatic invariants may be derived, and the will bounce back and forth 
across the line y = 0 with reflection in each of the arms of the channel. 

On the other hand, consider a particle starting its motion along the 
y-axis. By symmetry, the particle will maintain x = 0 at all times, and 
its motion will therefore be governed by the potential 

V = ( cr 2 2 /2ni) ( y - y Q ) 2 (76) 

which is a harmonic potential well. Presumably the orbits discussed in 

the section on ''intuitive treatment" will lie close to this class. If 

the corresponding motion in the x-direction were also harmonic, of the 

same frequency, we could get circular orbits in the averaged motion in 

the s = 0 plan© . In view of Sonne rup' s numerical results, this is only 
an d ©van that only , 

approximately true^in certain case®, while in general the averaged motion 
in the y-direction and in the z-direction have different properties. 

(4) Generalization 

Assertion ? In general, if to the field of eq. (l) we add a magnetic field 
independent of z and with no component in the z direction 



then the motion of a charged particle in the combined field can be reduced to 
motion in a potential in the (x, y) plane, and the field lines in this 
plane follow the equip© ten tial lines. 
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Proof : Clearly we can choose z to be one of the Euler potentials here 

z - P 

Let o/ = (x, y) be the Euler potential conjugate to this. The 

Hamiltonian then has the form 

H = (l/2o) | P x 2 + P y 2 + ( P z - (e/c) oC( x , y) ) 2 J* (78) 

Since z is absent, p z is a constant and our potential is 

V » (l/2m) ( p z - (e/c) <<(x, y) ) 2 (79) 

Magnetic field lines in the (x,y) plane stay in this plane by virtue of 
equations (77) and are characterized by constant values of ©4(x, y) . 

By equation (79) the potential V is also constant along such lines, 
hence Q.E.D. 

It is interesting to compare this case to Stdrmer' s treatment of 
particle motion in a dipole field. There, too, the problem reduces to 
motion is a two-dimensional potential field (Stdrmer: The Polar Aurora, 
Oxford 1955 j A.J.Dragt, Reviews of Geophysics 255, *65) , but the 
equipo tentials (except for V = 0 ) do not follow field lines. This may 
be traced to the fact that in Stdrmer's case the basic geometry is 
axisymmetric rather than independent of one of the cartesian coordinates. 
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X - t 


Null Point 


If to the main field of eq. (l) we add a field in the x-direction increasing 
linearly with y (i.e. a field resembling the main field in its form but 
orthogonal to it) 

B* - (b/L) y £ (80) 

then we obtain an x-type null at the origin (actually, a line of null points 
stretching along the z-axis). W© will assume that b << Bq , so that 
near the null point (80) represents a small perturbation. Since (80) is a 
field of the general type of (77) « it may be handled in the same manner as 
indicated there s 

B e = (b/L) y (vy x V z) = (b/2L) (Vy 2 X V z) 


V X (b/2L) y 2 z 


Xat us define 


0"^ S3 (]“ SS @ Bq/2Lc 


3b/ 2Lc 


B * (l/2m) £ p^ 2 + p y 2 + (p g + ^ x 2 - <T 2 y 2 ) j* 

As in (78), p is again a constant P of the motion. As in (7l)» the 

% 

Hamiltonian may be rearranged to the form 


2m ^ 


V * (* + 


* 2)2 } + t W * <».- «■***)*} 


"/2m - ( (r ± <rj- m) x y S 


(84) 
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As in (7l) » this may be viewed as representing two one-dimensional 

oscillators — anharmonic ones, of the "meandering” type — coupled by 

2 2 

a nonlinear term proportional to xy . This can be treated by pertur- 
bation methods in a maimer similar to what was indicated there. 


Alternatively, we can view ( 83 ) as representing motion in a two-dimensional 
potential 

T - k( p * +<r i * 2 - ^y 2 ) 2 (85 > 

Simultaneous inversion of the signs of (T^ and <T^ does not affect 

the Hamiltonian (although it also inverts the sign of P appropriate 

z 

to any given orbit), so the question is essentially whether the signs of 
these two coefficients are equal or not. Suppose that both are positive: 

<f 1 = 1 / S* (T 2 = 1 / S 2 2 (86) 


Then the equipo tentials for a given value of V will be two pairs of 
hyperbolas ? 9 




2 m V 


(87) 


As V becomes smaller the hyperbolas converge towards the pair 

( y 2 / S 2 2 ) - ( xf S 2 ) = P z (88) 

which corresponds to V = 0 ; these two hyperbolas are the bottoms of two 

potential valleys and are separated by a pass at which V = P /2m , which 
jalaoy z 

is^the potential on the two asymptotic straight lines . A sketch of the 
t for negative P j 

appropriate topographyCis given on the next page. If P_ is positive the 

z 

valleys of V *» 0 appear in the other family of hyperbolas, for only then 


can V = 0 be met with x = 0 



If QT^ and CT^ have opposite signs, we obtain a set of nested ellipses 
centered around an O-type null point. A single null point is not likely 
to be of this type, but if a chain of them is given (as in the next section) 
then O-type and X-type points alternate. 

C h ain of Null Points 


possibly 

A practical case, ^ resembling the tail of the earth’s field more than any of 
the preceding models, may be a chain of alternating neutral points of X-type 
and O-type, along the surface x = 0 . This was suggested in theoretical 
work of Coppi et al. (Phys. Rev. Letters 16, 1207, 27 June ’66) and by 
analyzing observations by Schindler and Ness (J. Geophys . Res. XL* 91* ' 72 ). 
The latter authors also provide a model for the field, but we will 
represent the situation here by a somewhat simpler model than theirs. Such 
a model is 

B = (Bq/l) x y + b x sin coy (89) 
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This may be written 

B = V x ( £o 008 00 y " (Bq/ 2L) x 2 ) £ 


leading to the Hamiltonian 


- ~ V x ( <T^ x 2 + (TgCOS Uy) z 

(90) 


H 




( p , + *i * 2 


+ <T 2 cos <Jy 



(91) 


This as before can be resolved into two one-dimensional coupled motions, 

^while/ 

of which one is oscillatory ~ T ' the other may or may not have a steady 
growth , in addition to its oscillation. 

Alternatively, we may view H as representing motion in a two-dimensional 
potential 

V = ~ (P z + 0"i x 2 + <T 2 cos Wy) 2 (92) 

This potential is non-negative and its gradient 


VY = T 0 (2^xx - sin Wy y ) . (93) 

vanishes if 

x = 0 ; oiy = iJTN (N = 0, 1, 2, . . . ) (94) 

The lines of constant V are sketched below: 
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Since for small y 

<r 2 COS coy £ <r 2 - ■£• < 3 " 2 CO 2 y 2 (95) 

we get hyperbola-like levels near the origin if 3"^ and ( are of the 
same sign, i.e. the origin and all null points with even N are of X-type. 

At these points 

* ■ k (p z + ^ <*> 

while for odd N we get O-type nulls and 

V ■ k (P z - r 2 )2 <*> 

Let us assume <T^ and <7* 2 are both positive (if both are negative, the 
results derived below all hold with the sign of everywhere reversed; 
if their signs differ, the whole pattern shifts phase). 

On the line x = 0 we have 

v = h (K + ^ cos 63 y ) 2 (98) 

dm z d 

If |P I < (T„ , (98) will pass through zero and therefore the potential 
z d 

will vanish on closed channels of roughly elliptical shape, each surrounding 
a central peak and separated from similar channels by two passes, the whole 
chain of passes and channels being flanked by two rising slopes. 

If | P I > (T- , the situation depends on the sign of P .If P > C T 

Z d Z Z Id 

the expression being squared in (98) is ad ways positive, and the extra term 

i 

2 

x* in (92) can only increase it. Thus V never vanishes and is smallest 
at the O-type nulls, where (97) holds, and these form a series of pits separated 
by passes, again the whole chain flanked by rising slopes. If P < - CT1 , the 

Z d 

same expression is always negative. However, by (92) the expression inside V 
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reverses sign and becomes positive for large lx j . We then get two wavy 
valleys of V = 0 flanking a central ridge of alternating peaks and passes, 
and these valleys in their turn are flanked as before by two rising slopes. 


\ 

\ 



